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Abstract 

We present an exact numerical computation of the one-loop correction of 
the string tension for the Nielsen-Olesen vortex in the Abelian Higgs model. 
The computations proceeds via the computation of the Euclidean Green's 
function for the gauge, Higgs and Faddeev-Popov fields using mode functions, 
and taking the appropriate trace. Renormalization is an essential part of 
this computation, it is done by removing leading order contributions from 
the numerical results so as to make these finite, and to add the divergent 
parts back, after suitable regularization and renormalization. We encounter 
and solve some problems which are specific to gauge theories and topological 
solutions. The corrections to the energy are found to be sizeable, but still 
smaller than the classical energy as long as is smaller than unity. 
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1 Introduction 



Extended classical solutions that can be interpreted as strings exist in various 
realistic and semi- realistic models of particle physics [H [2]. Their possible 
role in cosmology has been evocated long ago [3], see [H [5] for reviews. 
In confrontation with the recent wealth of cosmological observations they 
present a very interesting and active, but still controversial field of research, 
see e.g. P E]. 

In the present investigation we consider the string that is made up by the 
vortex solution of the Abelian Higgs model in 3 + 1 dimensions, well-known 
from superconductivity |B] , and commonly denoted in particle physics as the 
Nielsen-Olesen vortex [5]- The 1 + 1 dimensional version of this solution 
represents an instanton solution, that has been widely considered in the 
context of baryon number violation. 

There have been various investigations of the one-loop correction to the 
string tension. The fermionic corrections to the energy of the Nielsen-Olesen 
vortex have been computed exactly in Ref. [TDI. Fermionic corrections to 
strings in 2 -(- 1 and 3 + 1 dimensions have also been considered by various 
other authors [HI [12], |13l HI]. Such calculations may be important in the 
context of the instability of the electroweak string [15] 

Here we consider the energy corrections arising from the gauge- Higgs and 
Faddeev-Popov sector. There are some previous investigations, using heat 
kernel techniques [161 lU]- Here we attempt an exact computation, using 
techniques that have been developed previously in Refs. [TH [191 120] and 
applied in various semiclassical computations, of one-loop energy corrections 
and one-loop prefactors to transition rates. The computations have three 
essential ingredients: the use of mode functions ("Jost functions", see e.g., 
[2T] ) in order to compute exact results, the use of perturbative subtractions, 
so as to make these results finite, and the computation of the subtracted 
parts using a covariant regularization and renormalization scheme. Similar 
approaches have been used recently by other authors, see e.g. [221 l23l [211 
125]. In a gauge theory as considered here one finds some complications due 
to the fact that there are cancellations of divergences between graphs of a 
different number of vertices. This will discussed in detail in the context of 
renormalization. 

Whether the energy corrections are small or big depends, in the present 
case, of the gauge coupling. The classical string tension is proportional to 
= m'^/g'^, where g is the gauge coupling, while the corrections are pro- 
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portional to multiplied by a function of ^ = rriH/mw. So the corrections 
are necessarily small relative to the classical string tension if g is sufficiently 
small. If ever the corrections are big then this signals the breakdown of the 
semiclassical method. In intermediate situations one may have recourse to a 
Hartree type approximation by including the back-reaction of the quantum 
fluctuations to the classical solution. The methods we use here are suitable 
for such investigations EZ] , or even for self-consistent calculations without 
classical solution [28] 

The text is organized as follows: In section [2] we present the model, the 
classical vortex solution and the classical string tension. In section|3]we relate 
in general the fluctuation operator to the one-loop correction to the string 
tension, and we explicitly derive the fluctuation operator. Its partial wave 
reduction is presented in section HI This is the basis for the actual numerical 
computation of the one-loop string tension, which is described in section [51 
Renormalization is discussed in some detail in section [61 In section [7] we give 
some details of the numerical implementation and present the results. We 
conclude with a summary in section [HI Some technical details are discussed 
in Appendices A-G. 

2 Basic relations 

The Abelian Higgs model in (3+1) dimensions is defined by the Lagrange 
density 

^ = -If.^F^" + - ^ {\<j>\' - v')\ (2.1) 

Here is a complex scalar field and 

= d^A,-d,A^, (2.2) 
= d^-igA^. (2.3) 

The particle spectrum consists of Higgs bosons of mass m\j = 2Af ^ and vector 
bosons of mass = g'^v'^. The model allows for vortex type solutions, 
representing strings with a magnetic fiux, the Nielsen- Olesen vortices [U [9l 
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. The cylindrically symmetric ansatz for this solution is given by ^ 

Af'\x,y,z) = '-^A{r) z = l,2, (2.4) 



<P''{.x,y,z) = vf{r)e 



iip{x) 



(2.5) 



where r = a/o;^ + and if is the polar angle. Furthermore Af = Aq = 0. In 
order to have a purely real Higgs field one performs a gauge transformation 



-lip 



Ai ^ Ai-^i^lg 
to obtain the instanton fields in the singular gauge 



gj.2 



[A{r) + 1] i = l,2 



(t){.x,y,z) = vf{r). 



(2.6) 
(2.7) 



(2.8) 
(2.9) 



With this ansatz the energy per unit length, or string tension a takes the 
form 



TTf^ / dr 



rm 



w L 



dA{r) 



dr 



df{r) 
dr 



2 f2 

+ 



[A{r) + 1]^ 



+ 



The magnetic fiux is given by 



M 



d^xB^ 



dx dyFi2 



(2.10) 



(2.11) 



Explicitly we find 

= J d^xiV^A'^^-V^Af^) = j d(f)rdr 



— 1 1 2tt 

—A'{r) = — M(0) - A{oo)] 

gr J g 

(2.12) 

For the case run = an exact solution to the variational equation is known 
[29j . for which the classical string tension takes the value ad = vrf^. We here 



We use Euclidean notation for the transverse components, so = = —Ai etc. 



3 



will consider the general case m^y 7^ tuh, for which the classical equations of 
motion 

^ - - "•»./'('■)} + 11 = (2,14) 

have to be solved numerically. 

Imposing the boundary conditions on the profile functions 



A{r) const ■ , A{r) —1 

/(r) const ■ r , /(r) 1 



(2.15) 



the magnetic flux is $m = 27r/(yf, the Dirac magnetic flux quantum, and the 
action is finite. 

Since we have to consider fluctuations around these solutions a good 
numerical accuracy for the profile functions /(r) and A{r) is required. As in 
previous publications [311 [32] we have the method of Bais and Primack [33] . 
The values for the classical string tension are given in Table [H 



3 Fluctuation operator and one-loop string 
tension 

The fluctuation operator is defined in general form as 

^ ^ 5^;(a;)57/>,(x')''^''='^^' ' ^^'^^ 

where '^i denotes the fluctuating fields and '^f the "classical" background 
field configuration; here these will be the vortex and the vacuum configu- 
rations. If the fields are expanded around the background configuration as 
'^i = '^f + 0j and if the Lagrange density is expanded accordingly, then the 
fluctuation operator is related to the second order Lagrange density via 

C}' = \^\M.,<\>, . (3.2) 
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In terms of the flTictuation operators Ai on the vortex and on the 
vacuum backgrounds, the effective action is defined as 

o^^l A^. A^n , \ ■ (3-3) 



"^^^-2"^\detM0 + i,^ 

As the background field is time-independent and also independent of z the 
fluctuation operators take the form 

M = -dl\ dl - , (3.4) 

where is a positive-definite operator describing the transversal fluctua- 
tions. As is well known the logarithm of the determinant can be written as 
the trace of the logarithm. One can do the trace over pq, the momentum 
associated with the time variable, by integrating over T J dpo /27r, where T 
is the lapse of time. One then obtains 

= [E^ - £;f ] , (3.5) 

where are square roots of the eigenvalues of the positive definite operator 

-dl + M^, (3.6) 
and likewise Ea^ are those of the analogous operator in the vacuum 

- dl + = -dl - Vi + (3.7) 

Here = diag(m^, . . . ,m^) is the diagonal mass squared operator for the 
various fiuctuations. 

So the effective action becomes equal to the difference between the zero 
point energies of the fluctuations around the vortex and in the vacuum, 
multiplied by —T. We can also do the trace over the variable p3 by integration 
over L J dp^/2'n. We then obtain 



S., = -^TLY: J [^/^.-^^^' 



(3.8) 



where /^^ are the eigenvalues of the operator M.-^ and /Xa ^ those of — V^^ -|- 
m^. In the same way the classical action becomes 

S,i = -TLa,i (3.9) 
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where is the classical string tension. So our goal reduces to computing 
the one-loop approximation to the string tension given by 



O"l-loop = O'cl + 0"fl , (3.10) 

where the fluctuation part of the string tension is given by 



E[ dk3 1 
/ 2^2 



,{0) 2 



(3.11) 



Of course all expressions are formal, the integrals do not exist before a suit- 
able regularization. Anyway we do not plan to compute any eigenvalues 
of the fluctuation operators but will reduce these expressions to traces over 
Euclidean Green' s functions, where renormalization will be done properly. 
However the formal identities will allow us to trace the way in which the 
counter terms in the original Lagrangean enter the final expressions that are 
going to be computed numerically. 

The fluctuation operator has been derived previously [21], in the con- 
text of quantum corrections to the Abelian instanton; we here recall this 
derivation. The gauge and Higgs fields are expanded as 

= A^i + a^, (3.12) 
4> = 4>'' + ^. (3.13) 

In the following we will drop the superscript cl, so the letters A'^ and (p will 
denote the background field, and and ip the quantum fluctuations. 

In order to eliminate the gauge degrees of freedom we introduce, as in 
Ref. [31j, the background gauge function 

^(a) = 9X - f (0V - ^V'*) • (3.14) 

We note that for the background field (9^v4^ = 0. In the Feynman background 
gauge we get the gauge-fixing Lagrange density 

-^{dXY - ja''ifdf.<P + (pd^f - V>*d,<P - <pd,if*) (3.15) 
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The associated Faddeev- Popov Lagrangean becomes 

Cj.p^\r^\-d'-g^')r^. (3.16) 

In terms of the real components tp = ipi + i(p2 and rj = {rji + ir)2) the second 
order Lagrange density now becomes 

{C + jCgf + jCfp)" = -a^^ {-D + g^(f)^) 

+ ^2{9A^'d^)^i + ^i{-gA^'d^)^2 (3.17) 
+ a''(2/^^(/.)(^i+ a''{2gd^<P)^2 

(-° - 5^') ^1 + (-□ - 5V') ^2 , 

where we have omitted the superscript from 0*^ and A'jj). We now specify 
the fluctuating fields (■0i, ■02, ■03, ^04, ^^s) as (a^, a^, (pi, </72, 7712), 





( ai \ 


V'2 




03 


^1 


04 


<yi'2 


Us / 





where we have used Euclidean notation for the transverse gauge field com- 
ponents. We furthermore write (see also above) 

M = -d^o+dl-M^ , (3.19) 

separating the trivial part from the one that is modified by the background 
field. With these preliminaries we obtain the following nonvanishing compo- 
nents of the overall fluctuation operator A4j~^: 
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L 

11 








M- 


L 

22 








L 

13 


= 2g^Ai(l) 






M- 


L 

14 


= 2^Vr 






L 

23 








M- 


L 

24 


= 2^V2> 






L 

33 


= -A^ + /A^' 


+ g^(t>^ + 




M- 


L 

34 


= -9^- 






L 

44 


= -A^ + /A^' 


+ A(302 




M- 


L 

43 


= gA^- 


■V 



(3.20) 



As discussed above it is understood that the contribution of the Faddeev- 
Popov operator A^55 enters with a factor —2 into the definition of the one- 
loop string tension. The fluctuation operators for the vortex and vacuum 
background are now obtained by substituting the corresponding classical 
fields. The vacuum fluctuation operator becomes a diagonal matrix of Klein- 
Gordon operators with masses {■mw,'>Tiw,'>Tiw,'m'H,'mw)- It is convenient to 
introduce a potential V via 

M = M^ + V. (3.21) 
This potential will be specified below after partial wave decomposition. 

4 Partial wave decomposition 

The fluctuation operator Al"*" can be decomposed into partial waves with 
respect to the polar angle (p, and the string tension decomposes accordingly. 
We introduce the following partial wave decomposition for fields 



^ , , , f cos(^ \ e*"'^ . , , 
^ '\ simp J ^ > 

+ OC 



cos(p J 



/27r 

n=— oo 

+00 



n=— OO 
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After inserting these expressions into the Lagrange density and using the 
reahty conditions for the fields one finds that the following combinations are 
real relative to each other and make the fluctuation operators symmetric: 



Ff(r) 


= liUr) 


+ Cn{r)) , 






- c„(r)) , 


mr) 


= hn{r) , 






= ihn{r) , 






= 9n{r) . 





(4.2) 



Writing the partial fluctuation operators as 

= M^o + V , (4.3) 
the free operators M-'-q become diagonal matrices with elements 

(f Id nf 2 



where (rij) = (n— 1, n + 1, n, n, n) and (mj) = (m^, m^, m^, m^, mw)- The 
potential V takes the elements 







V" — 

* 12 — 






V 14 — 


^ 22 




* 23 ~ 


^24 




V" — 
*^33 — 




= -2 — 




^55 




V,5 = 



= 



= V2mwf 



A + 1 



+ 



m 



4^ 



) (/^ - 1) 



= 0. 



(4.5) 



Choosing the dimensionless variable m^r one realizes that the fluctuation 
operator only depends on the ratio mH/'mw up to an overall factor which 
cancels in the ratio with the free operator. 
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5 Computation of the effective string tension 



The method for computing the effective string tension used here is based 
on the Euchdean Green's function of the fluctuation operator. This Green's 
function is deflned by 

{u'' + kl + M^) g{x±,x'^,k3,p) = (5.1) 

and similarly for the operator A4^. It contains the information on the eigen- 
values of the fluctuation operator M.-^ via 



/ 



d^x^Tr g{x\ x\ h, u) = Y: o - (5.2) 



, >^l + kl + iy'- 
We define a function F{k3, u) as 

F{ks, u)= J (fx^ Tr {g{x^, x^, v) - Q^'ix^, x^, z/)) . (5.3) 



We then find 



(5.4) 



This expression is still to be integrated over ks and is by itself already linearly 
divergent. The sum over eigenvalues becomes an integral J d^k-^ and the 
difference of the energies behaves asymptotically as l/lfc-*-!. So regularization 
is required. This will be dicussed in the next section. Assuming that it has 
been achieved, we can do the v and k^ integrations at once, using the fact 
that F(/c3, v) only depends on p = ^kl + 1/2, i.e., F{k^,u) = F(0,p). We 
then obtain for the one-loop string tension 



After these more formal considerations we will present the way in which 
we actually compute F{k^,v). We first use the partial wave decomposition 
to write 



+00 



F{k^,v) = F{p)^ K{p): (5.6) 
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where 

F„(rt^/*.lV(G„(.nrt-G:(.nrt). (57) 

and where the partial wave Green functions are defined by 

(M„+p2)G„(r,r',p) = 1^5(r - r') . (5.8) 

For M° the Green function is simply a diagonal matrix with elements 

G° ,,(r,r',p) = 4^(«:,r<)K„^(«:,r>) , (5.9) 



where = y + p'^. For the Green function of the operator M„ the matrix 
elements similarly become 

G„ ,,{r,r',p) = f:-{r<)f:;{r>) (5.10) 

where the functions f^^ form a fundamental system of linearly independent 
solutions of 05.81) . regular as r — >■ for the minus sign and as r — oo for the 
plus sign. The correct normalization is obtained by imposing the boundary 
conditions 

f:tir) ^ 5tK^Mr). (5.11) 

as r ^ oo. Actually we have solved numerically the differential equations for 
the functions /i"^ defined by 

f:t = B^^{K,Tm^ + Kf{r)) (5.12) 

where i?+ = Kn and B~, = /„ are the appropriate Bessel functions, and 
with the boundary conditions /i" j — > as r ^ oo. In this way one keeps 
track of the free contribution oc 6T and 



Tr [G„(r,r,z/) - G°(r,r,z/) 
= W + /i^t W + Knr)h:nr)] InA>^^r)K^A>^,r) , (5.13) 



to be inserted into (15. 6p . 
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6 Renormalization 



Our numerical results are computed in such a way that we first compute 
the contributions of the subsystems for fixed n, then sum over n and finally 
integrate over p. The subtractions necessary to make the summation over n 
and the integration over p convergent are done in the partial waves. There 
are two ways of doing this: either one plainly subtracts in the partial waves 
the first order contributions TiGn^ = — TrG° V^jG^ and the second order con- 
tributions TiG'n^ = TiGlVijCiVjiGl, or one does the subtractions already 
in the functions /i"^ as described in detail in Refs. [20]. For the Faddeev- 
Popov sector we have used both methods and the results agree very well. In 
the gauge-Higgs sector the singularity of the external gauge field at r = 
leads to specific difficulties, to be discussed below. These are easier to handle 
with the first method, and so we decided to work with the plain subtractions 
throughout. 

The essential problems occuring here are twofold: (i) the external gauge 
field is not square integrable, due to its singularity at r = and (ii) there 
are cancellations between graphs of different order in the external vertices, 
the well-known cancellation of the quadratic divergence in the vacuum polar- 
ization involves seagull terms with one vertex and graphs with two vertices. 
Such cancellations occur in higher orders as well, so one has to be careful, 
some contributions with two external vertices do not have to be subtracted, 
because their divergences are cancelled in higher order. 

In the following we will treat the various contributions of first and second 
order term by term, comparing the contributions which are subtracted with 
the corresponding Feynman graphs. The subtracted contributions are then 
added back in covariant form. In this way we preserve covariance, which 
would be violated if we would introduce noncovariant cutoffs (e.g. in p ). 

The Feynman rules are formulated in the vacuum sector, with a gauge 
fixing analogous to the one for the vortex sector: we define the background 
gauge function 

J^,iA) = dX - f [(^ + V'OV - + (6-1) 
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and the gauge fixing Lagrangean 

= ~{^^,a>'y + ga^{f2d^,^i + {v + ifi)d^ip2) (6.2) 

2 I 2 I y 2\2 2 

-(y^^ + C/ W</?1 + y^i) <^2 • 

Then the Lagrangean, including the gauge fixing and Faddeev-Popov terms, 
takes the form 

+ ^d^^2d^'^2 - \m^w^\ (6.3) 
+g'^a^a^'v(pi + -g'^a^a'' [ipl + ip^) 

-g^vfifl - ^'^Ifl 

i=l 

The vertices for the Feynman rules are presented in Figs. [1] and [21 

We denote as Feynman graphs the vacuum graphs with the external 
sources which are provided by the classical Higgs and gauge fields: <fi{x) = 
v[f{r) - 1], (p2 = 0, ttf, = for /i = 0, 1 and = eikXk{A{r) + l)/gr'^ for 
i = 3,4. Vacuum graphs with external fields ip2 will not be displayed as they 
do not contribute. 

The vertex graphs originate from the expansion of the Green's function. 
Starting with 

(-A2 + A^o)G'(x,x') = -K(x)G(x,x') , (6.4) 
(-A2 + A^o)G'°(x,x') = 52(x_x') (6.5) 
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Figure 1: Vertices for Feynman graphs derived from the Lagrangean ( 16.3^ : 
gauge-Higgs sector. We have not included combinatorial factors for permu- 
tations of the external lines. 
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Figure 2: Vertices for Feynman graphs derived from the Lagrangean fl6.3p : 
Faddeev-Popov sector, without combinatorial factors for permutations of ex- 
ternal lines 



Figure 3: The one- vertex graph. The cross symbolizes the point where the 
first order Green' s function is contracted. 

we get the expansion 



G(x, x') = G°(x, x') - / (Px"G'>{x, x")r(x")G'°(x", x') + ... (6.6) 



which we recall mainly in order to keep track of the signs. The vertex graphs 
are obtained by contracting the Green's function according to Eq. (15. 7p in 
the partial waves, or, after summation, according to Eq. (15.31) . The first 
order and second order graphs are depicted in Figs. [3] and HJ respectively. 

6.1 Graphs with one vertex 

The graphs with one vertex all come with in a combination analogous to 
the one displayed in Fig. [5l except for the graphs with the seagull vertex. 
The external lines v(f)i and (pl combine as v'^{f{r) — 1) + v'^{f — 1)^/2 = 
f ^(/^ — l)/2. They then are in correspondence with one of the vertex graphs 




i 
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Figure 4: The two-vertex graph. The cross symbohzes the point where the 
second order Green' s function is contracted. 



with the external potential l^j(r). In order to illustrate the relation between 
the Feynman and vertex graphs we consider the graphs with a Higgs field 
loop. The corresponding vertex graph is the one with the vertex V44. We use 
the subscript 44 for identifying this contribution. 
The Feynman rules yield the contribution 



where = 2/e — 7 + ln47r and where we have used Xv"^ = rn?jj/2. The factor 
LT comes from the trivial integrations over time and along the string axis. 
The Feynman graph constitutes a correction to iS, where S is the action, so 
as an energy correction it receives a minus sign. 



Figure 5: Feynman graphs corresponding to the one- vertex graph with V4 
The contribution we subtract is given, in the partial waves, by 



= - / rdr r' dr'Vi (r)G° (r, r', (r', r, k„) . (6.8) 
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The corresponding vertex graph is of the type shown in Fig. [21 As announced 
above, we do not include the {A{r) + l)^/r^ term and denote the restricted 
potential by V^. We also leave out the index n as this potential does not 
depend on it. The r' integration can be carried out and one finds 

d f°° 

r drVu{r)Gl{r, r, kh] 



n,44 



dml JO 



(6.9) 



This can be summed up to yield 

d 



dmjj 
1 d 



271 



rdrVi{r)G\^,^, Kh 
d^k 1 



2pdp J (27r)2 + p2 + 

This is discussed in some more detail in Appendix 
dpp^/4:7i leads, after integration by parts, to 



d^xV^^r) (6.10) 
IBI Integration over 



Act 



fi,44 



p^dp 



An 
dudksd^k 



pdp d^k 



An (27r)2p+p2 + ^2^ 



d xVa 



44 I 



k2 + A;2 + Z/2 



44 



(6.11) 



{2nY + ftg + z/" + mj^ 

If we use dimensional regularization applied to the Euclidean four momentum 
integration dudk^d'^k = d'^k we obtain 



Aafl 



m 



44 



H 



32n^ 



Le + 1 - In 



m 



H 



m 



H 



647r2 



In 



m 



H 



d^V^^r) 



d'x{f{r) 



(6.12) 



in agreement with the result from the Feynman graph. In this case it was 
possible to check in detail the relation between the subtracted part, the vertex 
graph and the Feynman graphs. 

Obviously the procedure is analogous for the other graphs with one vertex. 
We obtain, both by using the Feynman graphs of Fig. [6] as also the vertex 
graph: 



A(T- 



fi,33 



H 



327r2 



+ 1 - In 



/i2 



d^xV. 



33 \ 



(6.13) 



— {mjj + 2m 



m 



w J 



B_ 

647r2 



+ 1 - In 



w 



d^x{f{r)-V) 
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Figure 6: Feynman graphs corresponding to the one- vertex graph with Vgf . 

For the contribution of Vn and V12 in the sector of angular momentum n one 
notes that they involve the Green's functions G^^i- In the sum the shifts 




Figure 7: Feynman graphs corresponding to the one-vertex graphs with Vn 
and V22. 

compensate each other and the result is, both from the Feynman graphs in 
Fig. [7] as from the vertex graphs, 



m 



1+22 



w 



327r2 



+ 1 - In ■ 



w 



m 



w 



167r2 



+ 1 - In 



m 



w 



d^x[Vu{r) + V22{r)] 
(Px{f{r) - 1) . (6.14) 



These are cancelled exactly by the Faddeev-Popov graphs of Fig. [81 




Figure 8: Feynman graphs corresponding to the one- vertex graph with V55. 
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6.2 Graphs with two vertices: Va x Va 

We have subtracted in the partial waves the second order contributions to 
^n\p) of the type 

Fn,ijAp) = j dr'r' j drr j dr" r"Vij{r)Gl{r,r\ Kj)Vji{r') 

G°(r',r",«:,)G°(r",r,«:,), (6.15) 

as represented graphically in Fig. HI The r" integration can be done, see 
Appendix [Cl with the result 

Fn,,Ap) = J dr'r' J r l-,^(r)G°(r, r', «:,)\/?^(r')G°(r', r, . 

(6.16) 

We recall that Kj = p"^ + rrij . The one-to-one relation with Feynman graphs 




Figure 9: Feynman graphs corresponding to the vertex graph with x 

is most obvious for the diagonal parts where i = j, and we will treat these 
first. For V33 and V44 we again discard the gauge field parts V^^ = V^^ = 
{A{r) + l)/r^, restricting them to the Higgs field contributions and 
This will be justified later. As now /tj = we can apply the derivative with 
respect to to both Green's functions and compensate this double action 
by a factor 1/2. Furthermore, we can replace this derivative by the derivative 
with respect to dp^ = 2pdp, applied to the whole graph. Then 




Figure 10: Feynman graphs corresponding to the vertex graph with Vgf x V3 
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P Fn,ii,ii{p) 



An 

1 f p^'dp 1 d 

2 J 47r 2pdp 



dr'r' / drrV^(r)Gl{r,r',K,)V^(r')G'i(r',r,K,) 



= -J^J^''''J drrV^{r)G'^{r,r',K,)V^{r')G'^{r',r,K,) . (6.17) 

As the potentials under consideration do not depend on n we may sum over 
n to obtain 

- j^-^ j ^^""j ^Vy,,(r)G°(x,x',«,)y,,(r')G°(x',x,«0 , (6.18) 

and this is just the conventional Feynman graph with external sources. As 
in the case with one vertex the integral over pdp/2n can be rewritten as an 
integral over dk3du/(27rY, and the Green's function G^ in two dimensions 
involves the integral over dk"^, so altogether we have an integration over 
d'^k /{2t:Y, times a factor 1/4. In momentum space we obtain 



_1 ! v.( \ ' I 



(27r)4[A;2 + m2][(A; + g)2 + m; 
where q = (0, 0, q) has only the transversal components. 



(6.19) 






Figure 11: Feynman graphs corresponding to the vertex graphs with Vu x Vu 
and V22 X V22 

The Euclidean Feynman integral is logarithmically divergent. In dimen- 
sional regularization it is given by 



J (2^ 



[F + m2][(A; + ?)2 + m2] 



(6.20) 



— a;)q^ + umj + (1 — a;)m^ 
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Figure 12: Feynman graphs corresponding to the Faddeev-Popov ghost ver- 
tex graphs with V55 x V55. 



This is to be folded with the Fourier transform of — 1: 

poo 

[f - l]lq) = 27r / drr{f\r) - l)Uqr) ^ xp-M , (6-21) 
Altogether we find 

Aafl,44,44 = I T^Al,- I d^^^ § \xp-ii<l)f ■ 



2567r2 J (27r)2 { " Jo 1^ 

(6.22) 

The corresponding Feynman graphs are displayed in Fig. O The external 
legs combine as 

18AV [{fir) - l)(/(rO - 1) + (/(r) - l)2(/(r') - l)/2 
+ - l)(/(r') - 1)V2 + (/(r) - l)^(/(r') - 1)74] 

= (9/2)AV(f(r)-l)(f (r')-l) (6.23) 

= (9/8)m^(/7r)-l)(/7r')-l). 



Using the MS* scheme the finite contribution is given by 
^^^MA^ = 1^J T^J d^^^- ^\xp^i{q)? ■ (6.24) 



2567r2 7 (27r)2 /i2 

The integral over the logarithm is discussed in Appendix O The Fourier 
transform of the potential, which actually is a Fourier-Bessel transform, is 
discussed in Appendix [F] The analogous diagonal contributions, displayed 
in Figs. [To], [11] and [12] are given by 

Aafl,33,33 = ^7-^ / 77^ / rf^ln \xp-i{q)\ , 



647r2 J (27r)2 

(6.25) 
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(6.26) 

and 

f (Pq r\ uj{l - uj)q^ + mfj 2 
Aa,,55,55 = -^ y ci^ln Ixp.^m , (6.27) 

repectively. The latter two contributions, those of the transverse gauge field 
and the Faddeev-Popov fluctuations, cancel each other, as for the flrst order 
contributions. 

6.3 Graphs with two verticesrVia x ^31 and V23 x V32 

The vertex graphs with x V31 and V23 x V32 correspond to the Feynman 
graph of Fig. [121 The Feynman integral is given by 




(6.28) 



p 



Figure 13: Feynman graph coresponding to the vertex graphs with x V31 
and 1^23 X ^^32- 

The external legs are given by f (/ — 1). Denoting the Fourier transform 

by 

POO 

Hf - l)]~(q) = vXf-M = ^-nv / drr{f{r) - l)Uqr) (6.29) 

Jo 

we obtain for the Feynman graph including the external legs 

'^^8^ J (2^^ l^^-^^^^l |L,-y rfcln 1 . (6.30) 
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The infinite part corresponds to the wave function renormalization of the 
Higgs field. In the MS scheme the finite correction to the string tension is 
given by 



A(Tfl,(i2)3(i2)3 = J \Xf-i{q)\ J duln 



w 



(6.31) 



6.4 Graphs with two vertices: Vu x V41 and V24 x Va 
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These nondiagonal terms correspond to the Feynman graphs of Fig. [TH The 
external Higgs field legs combine as 




Figure 14: The Feynman graphs corresponding to the vertex graphs with 
Vi4 X V41 and V24 x V42. 



1 + (fir) - 1) + (/(rO - 1) + (fir) - l)(/(r') - 1) = /(r)/(r') . (6.32) 
The Feynman integral is given by 

-^^-i^r^-j/o 7 r ^ 

The external legs combine to A^^, which restricts to a two-dimensional vector 
A^v f{r). We write the Fourier transform as 

[Aivf]~=27r-e,,q, / dr{A{r) + l)f{r)J,{qr) = -e,,q,XAf{q) ■ (6.34) 
9 Jo 9 

The Feynman graph with external legs then is given by 

■ rnn^w f fw^fr f\ 1 ^(1 - + ^^w + (1 - ^)^H \ 

^^^8^7 (2^ I r^-io 1 

(6.35) 
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Figure 15: Vacuum polarization Feynman graphs. They are equivalent to 
the vertex graphs with V34 x V43 and V43 x V34, and to the one- vertex graphs 
with and I44, respectively. 



The infinite term corresponds to the coupling constant renormalization in 
the term g'^\(f)\'^Afj_A'^/2. For the finite correction to the string tension we 
obtain 

'^O'fl, (12)4(12)4 = (6.36) 



Jo 



2 



87r2 J (27r)2 ''^^^'^^^ ^ 

6.5 Graphs with two vertices: V34 x V43 

The graphs containing V34 x V43 were the hardest stumbling block for this 
computation. As is well known, the quadratic divergence of the second order 
Feynman graph of Fig. [T5] cancels with the one of the first order seagull 
diagrams displayed in the same figure. The problem we have here is the 
fact that the quadratic divergence is proportional to A^A^ and this has a 
nonintegrable singularity at r = 0. While the ultraviolet divergence appears 
only after summation over n and integration over p, the radial integrations 
already appear in the partial waves. On the other hand, all we have to ensure 
is the fact that the final ultraviolet divergence is a wave function renormal- 
ization and does not contain terms proportional to A^A^. Furthermore our 
subtraction must render the total result finite. The solution of this problem 
was found essentially by trial and error and amounts to subtracting in the 
partial waves a second order term 



r dr r' dr'Vuir)Gl{r,r' , Kw) 



2 



rdr I r'dr'2n^^i^i^^G'i{r,r',Kw) 



^ [G'„(r ,r,/€H)J . (6.37) 
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and analogously for V43 x V34 or n^y <-> Kh- This contains already the sub- 
traction of first order divergences that would arise from the gauge field parts 
of V33 and V44. So the numerical results get finite without subtraction of those 
first order terms. However, this subtraction does not correspond exactly to 
the subtraction of the Feynman graphs of Fig. [151 it only takes care in the 
correct way of the divergent terms. The finite terms will be different. So 
when we add back the subtracted term after covariant regularization, we 
have to compute it in the same way as we did the subtraction, not by eval- 
uating the Feynman graphs with the external gauge field legs. This we will 
describe now. 

The Feynman integral corresponding to the second order diagram is given, 
in our gauge, by 



9 



d'^p 



(27r)4 (p2 - m]j + ie){{j) - qf - + ie) 



(6.38) 



where is the external momentum. The graph evaluates to 



47r2 



duj 



2Q2 



In 



2? 



+ q^q^uj <^ - In 



2! 



(6.39) 



with = ujm'^ + {l — 00)^2^^—00(1 — Lj)q'^. Our gauge potential is transverse, 
so in the following we do not have to keep track of the qf^q^ terms. The 
external fields as we have used them in the subtraction can be written as 



A{r) + 1 A{r') - A{r) 



gr gr 
A{r) + lA{r') 



(Air) 



lY 1 



gr 



9r' 



gr 



9r' 



(6.40) 



Taking the Fourier transform and setting q = q' because of momentum con- 
servation we have 



;eikqktjiqi{2'K) 



drJi{qr){A{r) + l) / dr' Ji{qr'){A{r') + 1) 



drJi{qr){A{r) + lY / dr'Ji{qr' 



(6.41) 
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Contracting with g^^, we have 



drJi{qr){A{r) + l) J dr'Ji{qr'){A{r') + 1] 



- j drJi{qr){A{r) + lf / dr' Ji{qr') 



(6.42) 



At first we consider the divergent part of the Feynman graph. It is given by 



9' 



IGtt 



1 2 

3^ 



(6.43) 



This is to be multiphed with the external leg factor given above and to be 
integrated with d'^q/{2TT)^. The part independent of q^ vanishes upon this 
integration: 



q dq 



drJi{qr){A{r) + l) J dr'Ji{qr'){A{r') + 1) 



- J drJi{qr){A{r) + lf J dr' Ji{qr') 

dr{A{r) + 1) J dr'{A{r') + 1) 
-J dr{A{r) + lf j dr'^ 



5{r - r') = , 



where we have used 



/ 



1 



qdqJn{qr)Jn{qr') = -S{r - r') 



(6.44) 



(6.45) 



This cancellation implies, as it should, the absence of a divergent term pro- 
portional to Af^A^^ in our subtraction. The remaining divergent term is pro- 
portional to g^. We note that now = — q^. Furthermore, J drJi{qr) = 1/q, 
J dqq'^Ji{qr) = and qJi{qr) = —dJo{qr)/dr. Including all factors we find 

^ ^ drJ^(qr)(A(r) + l) J dr' J^(qr')(A(r') + 1) 



167r2 3 



-2tt / q^dq 



~j drJi{qr){A{r) + lf 

. U I r , f, dJojqr) 



{A{r) + 1)J 



dr 



,d Jo{qr') 
dr' 



, Le 1, 
' 167r2 3 



27r 



dr 



[A'{r)f 



{A{r') + 1) 
(6.46) 
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This is just the kinetic term of the gauge field, multiphed with the wave 
function renormahzation factor {g^/3)L^/8n^. In the MS scheme this term 
is dropped. We are left with the finite part. Denoting 



27t J drJ,{qr){A{r) + 1) = XAiq) , 
2n [ drJ,{r){A{r) + l)^ = XA^{q), 



(6.47) 
(6.48) 



the finite part is given by 
1 r d^q 



3434 



duoQ^ 



ln^-1 



27r 



Xa{(i) XA^{q) 



(6.49) 

with Q2 = ujm?^ + (1 — (^)rn'^ + — uj)ci^. The integration over duj can be 
done analytically, we do not display the somewhat lengthy result. 



6.6 Graphs with two vertices: gauge field contribu- 
tions in V33 and V44. 

There are some contributions with two vertices which we have left out. De- 
noting the gauge part of V33 with V33 with V33(r) = {A{r) + l)'^/r'^ and anal- 
ogously for V44 we have left out second order the diagonal parts V33 x Vgf 
and V33 X VJj3,and likewise for the 44 (Higgs) channel. These diagrams with 
seagull type vertices combine with higher order graphs to give finite results. 
So they do not have to be subtracted. The corresponding Feynman graphs 
are displayed in Fig. [12] and [T71 If one looks at the Feynman graphs there 
are graphs with more than three vertices that are superficially divergent. In 
fact these combine in such a way that their divergences are cancelled. So 
with the subtractions described in the previous subsections we have done 
the necessary steps towards computing the finite result. Indeed with these 
subtractions the numerical results get finite, and we have already presented 
the finite expressions by which these graphs are to be replaced. 



7 Numerical results 

We have carried out the numerical programme as described in the previous 
sections, for values of ^ = mu/mw between 0.5 and 2. 
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Figure 16: Feynman graphs corresponding to the vertex graph with VJj| x 
and the three-vertex graphs that cancel their logarithmic divergences. 
The cancellation for the graphs related to I44 x is analogous, with the 
replacement 0i in the internal lines. 

The programme starts with computing a fundamental system of the mode 
functions h'^{r) for the 4x4 gauge-Higgs and the Faddeev-Popov sector. 
We have used 2000 grid points in r, up to rmax = 30, as we did already 
for the classical profiles /(r) and A[r). These are the basis for the partial 
wave Green's functions. This computation is identical to the one performed 
already in Refs. [311 [32]. 

We have computed the partial wave Green's function and the related 
integrals up to n = n = 35. Contributions of higher n were included by fitting 



Figure 17: Feynman graphs corresponding to the vertex graphs with V33 x 
1/33 and V^4 X Vl^ and the higher order graphs that cancel their logarithmic 
divergences. 
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the data between n — 5 and n using power fits An^^ + Bn~^ + Cn~^, and by 
appending the sum for n < n < oo on the basis of these fits. The perturbative 
subtractions were done in the partial waves, the sum over the unsubtracted 
and the subtracted partial wave contributions constitute the unsubtracted 
and subtracted functions F{p). These are displayed, for ,^ = 1 in Figs. [18] 
and [ini for the gauge- Higgs and the Faddeev- Popov sector, respectively. All 
the general features of these figures are similar for other values of ^. In these 
figures we also present the perturbative contributions of first and second 
order, summed separately. These can and have been used for cross-checks 
against semi-analytical results: 

The numerical sum over the first-order subtractions can be checked against 
the result obtained by summing the partial wave contributions analytically. 
One finds 



=1,4 



3ml, + m\j/2 



9k2 



+ 



3m^/2 
2k\ 



r dr{f'^{r) — 1) 



(7.1) 



The Faddeev-Popov contribution, which was computed separately, behaves 
(including the factor -2) as 



m 



w 



''W 



r dr[f'^{r) 



11 



(7.2) 



The integrals over r, here and below, are computed numerically using the 
classical profiles. 

For the second-order contributions we cannot do the sum over partial 
waves analytically, this involves the numerical Fourier transforms of various 
functions of the classical profiles; but we can easily check the asymptotic 
behaviour. It is given by 

E%,(p) = {mty + {ml, + mj,/2Y + 9m%/A) Jrdrifir)-!)' 



+8m^ / rdrf'{rY + Sm^ j rdrp{r)- — ^^y- — ^— 



+4 r dr 



[A'{r)f 



1 



(7.3) 



The first term collects the diagonal contributions that can be derived analyt- 
ically by summing up the partial waves. The three other contributions are 
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retrieved from the Feynman graphs by doing the transverse loop momentum 
integrations only. They collect the 1313 and 2323 contributions, the 1414 and 
2424 contributions, and the 3434 contribution, respectively. Of course these 
can also be computed and verified separately. The asymptotic behaviour of 
the Faddeev-Popov contribution is given by 

i^5555(p) = 2^/ rdrifir) - if . (7.4) 

These semianalytic results can be and have been checked against the numer- 
ical subtractions, thus verifying prefactors and signs. 




Figure 18: The integrand function F{p) defined in Eq. (15. 6p . for the gauge- 
Higgs sector: circles: the unsubtracted function; dashed line: one vertex 
contribution; dash-dotted line: two-vertex contribution; squares: subtracted 
function; straight line at small p: translation mode pole 2/p^; solid line at 
large p: asymptotic fit. 
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10 1 




Figure 19: The integrand function F{p) defined in Eq. fl5.6p . for the Faddeev- 
Popov sector. Symbols and hues as in Fig. [181 except that there is no 
translation mode pole at small p. 

After the subtractions the function -Fsub(p) behaves as and the integral 
over p'^dp can be done in order to obtain the subtracted part of the string 
tension. Finally we have to add back the subtracted terms in a covariantly 
regularized and renormalized form given in the previous section. 

As we have already mentioned, the integrands for the p integration are 
displayed in Figs. [18] and [191 ^^i the gauge-Higgs and the Faddeev-Popov 
sector, respectively. The figures show the unsubtracted functions, the first 
order and second order contributions and the subtracted functions. Note 
that the functions differ, at large p, by several orders of magnitude, so these 
subtractions are quite delicate. 

The gauge-Higgs sector displays a behaviour for small p, which is 
of course related to the 2 translation modes of the two-dimensional solution 
(the "instanton"). Here it merges into a continuum of transverse string os- 
cillations. The small-p pole remains of course after subtractions and causes 
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the gauge Higgs sector to be much more important than the Faddeev-Popov 
sector. Indeed, the finite parts obtained after integration with p^dp/ATi are 
much larger for the gauge-Higgs than for the Faddeev-Popov sector. Obvi- 
ously this important contribution of the translation mode is related to the 
transversal quantum oscillation of the string. 

At large p the subtracted integrand should behave as p~^. For the 
Faddeev-Popov sector this is realized in ideal form. For the gauge-Higgs 
sector there are deviations for p > 6, which remain even if higher partial 
waves are included. So they seem to be caused by some very small numer- 
ical deficiencies in the low partial waves Their origin is difficult to localize. 
As these contributions only appear above p ~ 6 we have done least-square 
fits of the form A/p^ + B/p^ based on the data points between p = 2 and 
p = 6. These fits then were used in order to append the integrals from p = 6 
to oo. The fits are displayed in Figs. [18] and [19] for the gauge-Higgs and 
Faddeev-Popov sectors, respectively. 

The integrals over the subtracted part of the function F{p), including the 
asymptotic tail based on the fits are given in Table [T] they are denoted as 

Now that we have computed the subtracted integrals we have to add 
back the regularized and renormalized divergent contributions. Those with 
one vertex are given, in unrenormalized form, by 

Aa« = {(L, + 1 - In ^)(1 + 2^^) - le^ne^ J dHf{r) - 1) . 

^ (7.5) 
One may choose the renormalization such as to omit these tadpole contribu- 
tions entirely, as one would do in the 1-1-1 dimensional theory, using normal 
ordering the field operators. Here we use the strict MS scheme and choose 
the renormalization scale as /i^ = m^. Therefore we have to add back 

A4:,) = -^{(l-ln!^)(l + 2a-^^^ln^^}|dMf(r)-l) (7.6) 

As we separately present this contribution in Table [I] the reader may easily 
change this contribution according to her/his own preferences. 

The separate results for the graphs with one vertex, ^<7^l, for the graphs 
with two vertices Acr^^ , the subtracted contribution AcTgub and the total 
one-loop contribution Aatot are listed in Table [1] and displayed in Fig. [20] . 
We should like to recall that the subtracted contribution still includes some 
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gauge field graphs with two vertices, which by themselves are divergent, but 
whose divergences are cancelled by higher order vertex graphs, as discussed 
in subsection 16.61 




Figure 20: The correction to the string tension as a function of ^ = mH/mw. 
Squares: contributions with one vertex, Ao"'-^''; triangles: contributions with 
two vertices Acr^^-'; diamonds: subtracted contribution Aarmsub, circles: total 
one-loop correction Ao"tot, asterisks: the classical string tension mutiplied by 
(yf^/vr. All string tensions are in units of m^. 



8 Summary 

We have presented here an exact, albeit numerical, computation of the one 
loop corrections to the string tension of the the Nielsen-Olesen vortex in 



33 



the 3 + 1 dimensional Abelian Higgs model, taking into account the fluctu- 
ations of gauge, the Higgs and the Faddeev-Popov fields. One of the main 
complications arouse from the fact, that the divergences of graphs with exter- 
nal gauge field loops require cancellations of Feynman graphs with different 
numbers of external vertices. This has required an extensive discussion of 
vertex graphs and vacuum FejTiman graphs. On the numerical side, we had 
to adapt a previously developed computation scheme to this new situation. 
We had to find and have found a way of implementing the necessary cancel- 
lations, which are relatively straightforward when done analytically, into the 
numerical procedure. 

The size of the corrections of course depends on the renormalization 
scheme and renormalization conditions. We here have adapted the MS 
scheme with TTly^ clS the renormalization scale. The corrections are sizeable, 
but small with respect to the classical string tension as long as the gauge 
coupling g is smaller than unity, which may be considered as a reasonable 
assumption. Unlike in the case of the fermionic corrections we do not have at 
our disposal an extra parameter, the Yukawa coupling, that could render the 
corrections important for heavy fermions. The parameters are fixed already 
at the classical level. 

Of course, in the present situation of cosmic string phenomenology a 
precise information on the string tensions and corrections to it cannot be 
considered as very important. Indeed we consider as our main result that we 
presented a method for computing such corrections in a gauge theory with 
all its technical complications. The situation may be different if these correc- 
tions are computed at finite temperature, as it may be realistic in primordial 
cosmology. Such computations, using similar techniques, have been done for 
bubble nucleation in Ref. [35,j. In the high-temperature approximation to 
the electroweak theory the corrections to the transition rate were found to 
be huge [Ml [37]. 

A further application of the methods presented here may be the inves- 
tigation of the role of quantum fluctuations for the electroweak string, in 
particular in the context of its stabilization [38] 
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Appendix 

A Partial wave decomposition of the free 
Green's function 

We shortly recall the partial wave decomposition of the free Green's function, 
a decomposition that is used repeatedly in the subtraction procedure. We 
here consider the Green's function in two dimensions, with an effective mass 
= m^+p^; the one in four dimensions is obtained by replacing = z/^ + /c| 
and by further integrations over u and k^. We further omit the subscript _L, 
replacing x_|_ ^ x and k_|_ k. The free Green's function Gq{:k,x.' , k) is 
defined as 



It can be computed readily, using formulas (9.1.18) and (11.4.44) of Ref. [39] 




with R = |x — x'|. Furthermore this may be expanded, using Eq. (4) in 
section 7.6.1 of Ref. [IQ] as 




(A.l) 





(A.2) 




(A.3) 
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B First order contributions to the partial waves 

Using the formula for the free Green's function derived in Appendix |A] the 
first order contribution to the Green's function in two dimensions is given by 

G'(^)(x,x',p) = - J d^x"Go{^,^",K)V{r")Go{^",^\K) 

= I Yl 5^ ^n(/tr>)J„(/tr<)e™('^-^") 

n=— oo n'=— oo 

xV(r")K„'«)4'«)e^'^'^^""^'^ (B.l) 
= / dr"r" J2 i^n(/«r>)/„(/tr<)e^"(^-*') 

n=—oo 

X\/(r")i^n«)/n«) 

where r> = max{r, r"}, r^^ = max{r', r"} and similarly for r< and r'^. The 
"potential" V{r) subsumes all vertex contributions appearing in first order. 
We define 

oo 

G«(x,x',p) = - V e-(<^-^')G«(r,r',p) (B.2) 
Ztt ^ — ' 

n=— oo 

and therefore get 

Gll\r,r,p) = - [ dr"r"Kl{Kr^)ll{Kr^)V{r") . (B.3) 



Using the integrals of Appendix [D] this can be integrated over r dr with the 
result 



dr r G^^\r,r,p) = — dr'r'V{r') 



oo 



^0 

z 



x—{-2n [ln+iiz)I^{z)Kl{z) + Kn+i{z)Kn{z)I;[z 



< — /. I I I # 1 -I I X I # I X I r\ \ /. \ —I— r\ 1 1 I X I r\ „ I X I # ^ 
+^ W^^ W - W^n } , (B.4) 



with z = kt' . Using the Wronskian relation 



In{z)Kn+l{z) + In+i{z)Kn{z) = ^ (B.5) 
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this can be rewritten as 



dr r Gll\r,r,p) = - I dr'r'V{r')— {-2nIn{z)Kn{z) 



+Z [Kn+l{z)In{z) - In+l{z)Kn{z)]} (B.6) 

Id 

- / drVVir') ——l4z)Kr,iz) . 



This is the contribution of one partial wave to the sum over partial waves, 
which is then to be integrated with respect to v and to yield the string 
tension in first order. However, the sum over partial waves and the integrals 
are divergent; the perturbative contribution has to be subtracted from the 
full, nonperturbative contribution of each partial wave in order to obtain 
covergent summations and integrations. We may check that formally the 
first order perturbative contributions can be summed up and integrated so 
as to reobtain the Feynman integral. The sum over partial waves of this 
expression yields 

1 d f°° 

= -^""Tpd-pi drVVir')Goi^,.,.) 

^ d f ^r, f d^k^ 1 ,^ , 

= -YpTp J ^ J (2^kiwT7 • ^ 

Going back to section [5] we see that we have integrate this with — J p^dp/in. 
We obtain 

-'- f j2 TA/i f P^dpd^k±_ Id 1 



(27r)3 2p dz/ + p2 + 



1 /" 2 TA/i IN [ pdpd'^k± 1 



, rfWrixl) / £-£1— (B.8) 

1 ,2 TA/i IN dvd^k 1 



2 J " "7 (27r)4 k2 + z/2 + m2 ' 

where in the last step we have substituted p^ = k^ + u"^. The result is indeed 
the first order vacuum graph. 
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C Second order contribution to the partial 
waves 



Using similar steps as in Appendix [B] and using again the integrals in Ap- 
pendix [D] and Wronskian relations, we find that the second order contribution 
to the partial waves is given by 




where r< = min(r', r") and r> = max(r',r"). The computation is quite 
lengthy, as we start with an expression where the values r, r' and r" appear 
in six different orderings. Here we have written all relevant indices, as the 
two propagators have, in general, different masses and angular momenta. We 

note that the differentiation with respect to Kj = y^p^ + can be written 

as a differentiation with respect to p is the combinations VijVji and VjiVij 
are combined. If applied to the sum over i and j, this results in a double 
counting that has to be compensated by a factor 1/2. So 




x/„^(Kjr<)ir„^(Kjr>)/„^,(fi;jr<)ir„^,(/€jr>) . (C.2) 

D Some integrals 

In our calculations we repeatedly need some integrals over Bessel functions 
which we write down here, without derivation. They can be obtained from 
Eq. (11.3.29) of Ref. [39] by taking the limit k — >■ / (see there). The formulas 
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are: 



jyrrlliKr) = ^ {-2n4+i(^)4(^) + z [ll{z) - l'^^,{z)] } , (D.l) 

f°° z 
j drrKlinr) = — {-2nK^+,{z)K^{z) - z [KI{z) - Kl^,{z)] } ,(D.2) 

/ 



z 

drrIn{Kr)Kn{Kr) = — {n[Kn+i{z)Iniz) - ^^(z) 4+1(2;)] 

+z [In{z)K^{z) + /„+i(z)ir„+i(z)]} , (D.3) 



with z — Kr. Note the hmits of integration, the third integral is an indefinite 
one. 



E Some useful sums 

We define the two-dimensional Green's function 6*2°^ (x, x', i/^) as the solution 
of the equation 

[-A2 + m^ + u^] Gjo) = (^'(x - x') . (E.l) 

We readily obtain 

(x,x',i.2) = / , . (E.2) 

^ ^ ' ' ' j (27r)2 k2 + m2 + i/2 ^ ' 

The integral may be done exphcitely with the result 

Gf\^,^y) = ^K,{^iR), (E.3) 

with 



i? = |x - x'l = + r'2 - 2rr' cos{ip - ip') . (E.4) 
Furthermore the Gegenbauer expansion of the modified Bessel function yields 

00 

(x,x',i.2) = -Y, e-(^-^')4(«:r<)K„(«:r>) . (E.5) 

n=— 00 

The limit x — > x' of the Green' s function does not exist. However, what we 
need is 

In{^r)Kr.{Kr) . (E.6) 



2k dK 

n=—oo 
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We note that 

The hmit i? — > exists and one obtains 

1 d °° 1 

n=— oo 

F Fourier-Bessel transforms 

For the evaluation of the perturbative contributions we need the Fourier- 
transform of the external sources. Let us consider at first a scalar field 
$(rx), independent of Euclidean time r and z. In our application (j) will be 
f{r±) — 1, P{r±) — 1 etc.. We have 

^iQo,qz,q±) = J d^xe~'P"'^{r±) 

POO PTT 

= {2nyS{po)S{p,) / drr / d(^$(r)e-*«^^"°^'^ (F.l) 

Jo J-TT 

poo 

= {2T:f5{pMPz) / rdr^r)Uq^r) . (F.2) 
Jo 

In the following we will omit the trivial factor (27r)^5(po)^(Pz) originating 
from the r and z integrations and write 

/•oo 

^{q) = X*(9) = 2% rdr^{r)Jo{qr) (F.3) 
Jo 

implying q q±. The relation between the Fourier-Bessel transformation 
and the inverse transformation implies the relation 

J qdqJi{qr)Ji{qr') = -5{r - r') (F.4) 

for arbitrary /. We deduce the Parseval equation 

/ (0'^*^^^'' = / j^rdr\^{r)\' = j dM^r)\' . 

(F.5) 
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which can be used as a numerical cross check. If one disregards V34 = V43 
and the gauge field parts of V33 and V44 these Fourier-Bessel transforms are 
unproblematic. We shortly discuss the various components. 

In the diagonal we have the Higgs field parts which are proportional to 
/^(r) — 1. This function goes to zero exponentially as r ^ 00, so the Fourier 
transform exists. It is found to decrease exponentially at large q. 

Considering the contribution of V13 = V23 = V2mwf' the function /'(r) 
decays exponentially as r — >■ 00 and goes to a constant as r ^ 0. So Vi^^q) 
exists. It is found to decrease as 1/g^. 

The function Vu{r) = \plm-wf{^A^ l)/r goes to a constant as r — > and 
decreases exponentially as r — 00. So again the Fourier transform exists, 
and again it decreases as 

All the Fourier transforms discussed up to now have to be folded with the 
Feynman integral kernels (see Appendix [G]) which behave as In at large 
q. One easily convinces oneself, that the integrals over the kernels times the 
squared Fourier transforms are well convergent. 

The Fourier transform of the gauge field 

A- = .„%^<^ (F.6) 

is problematic. The function is not square integrable, due to the singularity 
at r = its norm diverges logarithmically. So the Fourier transform is not 
square integrable either. Nevertheless it can be computed. The Fourier 
transform must be of the form 

= . (F.7) 

We have 

XA(g) = eijqjA\{q) 

= e,,e,fc J £x^^^e-'%Xk (F.8) 

r dr ^ ^ / e~"'"'°'^cos(p 

[°° Air) + 1 

= 27r / r dr Ji{qr) . 

Jo r 

The function xa(q') behaves as /5g as g ^ with 

POO 

(3 = Tx j r dr{A{r) + 1) , (F.9) 
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and like 27r/g as r — oo. This asymptotic behaviour makes the square 
norm of the Fourier transform logarithmically divergent, as a reflection of 
the logarithmic divergence of the norm of the field Aj-[r). The handling of 
the graphs with two external gauge field legs is discussed in some detail is 
subsection 16.51 



G The second order kernels 

The finite parts of the second order graphs all contain an integral 

I{ml,mlq^)= [ du \og[u;{l - u)q^ + uml + {1 - u;)ml] . (G.l) 
We generally assume > 0, ml > 0. If furthermore > then we define 



= — V — V V • ^ ^ ^ 

One easily realizes that uo^ > 1 and ct;_ < 0. With this definition 1 is given 
by 

T{ma'^,ml,q'^) = uj+\nuj+ + uj^\n(—uj-) + {1 — uj+)\n{uj+ — 1) 

+ (1 -cj_)ln(l -cu_) +lng2 -2 . (G.3) 

This expression is symmetric in and m^. 

If = the integral is given, for ^ ml by 

Ami ml 0) = _ 1 . (G.4) 

— mi 

Finally, if ml = ml we have 

X{ml, ml, 0) = In ml. (G.5) 

For small the integral behaves as Ci + C2g^ where the constant Ci is 
of course identical to X^(m^,m^,0). For large the integral behaves as 
Ing^ + const. + 0(l/g^). 

The correction to the gauge boson propagator, the term Acrfl,3434 contains 
an integral 

J = [ dujQ^ 



1 r Q2 



(G.6) 
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with = Lomjj + (1 — oj)m'^ + — 6i;)q^. Using the same conventions for 
a;_ and a;+ with = my,, rrib — mn- It is given by 

J — (f [(9ci;_ + 9a;+ — 18a;_a;+ — 6) ln(— u;_a;+ + a;_ + a;+ — 1) 
+ (9a;_ + 9a;+ - 18a;_w+ - 6) In 

+3w!(3w+-tu_)lni^^+ (G.7) 

+3cc;^ (3ti;_ — a;+) In — 

-3c;! + l^uj-uj+ - 6cu_ - 3a;^ - 6tu+ + 4] /18 . 
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A 1 


A 2 


AcTsub 


AcTtot 




0.5 


.244 


-.001 


-.731 


-.488 


.75742 


0.6 


.204 


.002 


-.773 


-.567 


.81306 


0.7 


.176 


.005 


-.795 


-.615 


.86441 


0.8 


.154 


.009 


-.811 


-.649 


.91232 


0.9 


.135 


.014 


-.825 


-.675 


.95737 


1.0 


.119 


.022 


-.837 


-.700 


1.0000 


1.1 


.105 


.030 


-.848 


-.713 


1.0405 


1.2 


.092 


.041 


-.859 


-.726 


1.0792 


1.3 


.080 


.054 


-.870 


-.736 


1.1163 


1.4 


.069 


.069 


-.881 


-.743 


1.1518 


1.5 


.058 


.086 


-.891 


-.747 


1.1860 


1.6 


.048 


.107 


-.902 


-.747 


1.2190 


1.7 


.038 


.130 


-.913 


-.745 


1.2509 


1.8 


.028 


.155 


-.923 


-.740 


1.2818 


1.9 


.018 


.184 


-.934 


-.731 


1.3116 


2.0 


.009 


.216 


-.944 


-.719 


1.3406 



Table 1: The correction to the string tension as a function of = itih/itiw- 
We present the finite part of the contributions with one vertex, ^cr^l-, the 
finite contributions from graphs with two vertices Augj/, the sum of higher 
order terms AcTgub, and the total one-loop correction Aatot- We also include 
the classical string tensions. All entries are in units of m^. 
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